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The pyrochlore lattice involves comer sharing tetrahedra and the resulting geometric frustration is believed 
to suppress any antiferromagnetic order for Mott insulators on this structure. There are nevertheless short-range 
correlations which could be vital near the Mott-Hubbard insulator-metal transition. We use a static auxiliary- 
field-based Monte Carlo to study this problem in real space on reasonably large lattices. The method reduces to 
unrestricted Hartree-Fock at zero temperature but captures the key magnetic fluctuations at finite temperature. 
Our results reveal that increasing interaction drives the non magnetic (semi) metal to a ‘spin disordered’ metal 
with small local moments, at some critical coupling, and then, through a small pseudogap window, to a large 
moment, gapped, Mott insulating phase at a larger coupling. The spin disordered metal has a finite residual 
resistivity which grows with interaction strength, diverging at the upper coupling. We present the resistivity, op¬ 
tical conductivity, and density of states across the metal-insulator transition and for varying temperature. These 
results set the stage for the more complex cases of Mott transition in the pyrochlore iridates and molybdates. 


I. INTRODUCTION 

The presence of geometric frustration disfavors long-range 
antiferromagnetic order and promotes a complex magnetic 
state in insulating magnets with short-range interactioiP^. 
Two complications arise in magnetic insulators close to a Mott 
insulator-metal transition (IMT): (i) the ‘virtual hopping’ of 
the electrons mediate long-range and multispin coupling, and 
(ii) the magnitude of the local moments weaken as the sys¬ 
tem is pushed towards the IMT. The first effect can lift the 
degeneracy of the short range model and promote an ordered 
state while the second tends to destroy the magnetic state al¬ 
together. The outcome of this interplay is lattice specificP® 
of relevance to several real life materials, and requires tools 
beyond those usually applied in frustrated magnetism. 

The possibilities in charge transport are also interesting. 
While the large moment, gapped, Mott phase is insulating, 
the strong suppression of the moment near the IMT, and pos¬ 
sible orientational randomness due to frustration, can generate 
a ‘bad metal’ state on the disordered magnetic background. 
Such a state can involve a pseudogap, an unusually large low 
temperature resistivitjED and, possibly, an anomalous Hall 
response if the moments organize in a non coplanar manne]®. 

The pyrochlores are a fascinating structur^to explore these 
effects. In the deep Mott state on a pyrochlore lattice, where 
one expects only nearest neighbor antiferromagnetic coupling, 
the effective model can be written as a sum, over the tet rahe- 
dra, of squares of the total moment in each tetrahedrorP®! 
The minimum of this is infinitely degenerate since the four 
spins at the vertices of each tetrahedron just need to satisfy a 
zero vector sum. The appearance of longer range couplings as 
the electron-electron Hubbard repulsion reduces (or the band¬ 
width increases) can, potentially, lift the degeneracy and pro¬ 
mote some ordered state. The transport and spectral character 
between the spin liquid Mott insulator and the band semimetal 
is also not known. 

There are, remarkably, experimental realisations of Mott 
physics on the pyrochlore structure as observed in the rare- 
earth molybdatesi^Hni and iridate^^^^EI] These, however, in¬ 


volve 4d or 5d electrons and require Hund’s coupling or spin- 
orbit interaction, in addition to the Hubbard interaction, for 
a successful description. To set the stage for these more re¬ 
alistic but complex models this paper focuses on the role of 
just the Hubbard interaction on the half filled pyrochlore lat¬ 
tice. We address two broad questions: (i) What is the nature 
of the magnetic state as one moves towards weaker interaction 
in the Mott insulator, and in particular what are the magnetic 
correlations near the IMT? (ii) What is the impact of these 
magnetic correlations on electron physics (resistivity, optics, 
spectral features) near the IMT? 

Our main results, obtained by using a real-space static 
auxiliary-field-based method, are the following: (i) Increasing 
interaction in the ground state leads successively to a small- 
moment metal, then a narrow insulating pseudogap window, 
and finally to the gapped Mott insulator, (ii) The resistivity 
and the low energy density of states have a strongly non mono¬ 
tonic temperature dependence near the metal-insulator transi¬ 
tion (MIT), (iii) The finite moment phases, from near the MIT 
to the deep Mott regime, are all disordered. Well in the Mott 
phase they display what seems (within the limits of our sys¬ 
tem size) to be power law spatial correlations that survive to a 
small finite temperature. 

The paper is organized as follows. In the next section we 
provide a summary of existing results on the pyrochlore struc¬ 
ture - in particular the Heisenberg and Hubbard models. This 
is followed by a description of our method. The next section 
presents results on the phase diagram, density of states, trans¬ 
port, and optics, across the MIT. The final section expands on 
a few issues like effective models for the magnetism, compu¬ 
tational issues, and the magnetic structure factor. 


II. PREVIOUS WORK 

The pyrochlore Heisenberg antiferromagnet model, which 
describes the deep Mott limit, has been well studied, the Hub¬ 
bard model much less so. We summarize the key results from 
these and then move to our method and results. 






A. Pyrochlore Heisenberg model 


C. Hubbard with additional Interactions 


While the Hubbard model at half-filling and large interac¬ 
tion maps on to the S' = 1/2 Heisenberg antiferromagnet, we 
provide a more general discussion of this model below. 

The nearest-neighbor classical Heisenberg antiferromag¬ 
net, H = J > 0, on the pyrochlore 

lattice doe s not show any magnetic order down to zero 
temperaturJi®] (j’ = q). The model can be written as 
H = (J/2) J2a Pa’ upto a constant, where = Z]i=i 
where the are spins on the tetrahedron a. The ground state 
needs to satisfy the constraint Pq, = 0 for all a. This re¬ 
sults in an infinitely degenerate manifold of possible configu¬ 
rations {Si}. It has been established that there are no internal 
energy barriers between these degenerate minima at T = 0, 
and no free-energy barriers at finite T, precluding a freez¬ 
ing transition eithei^^. Since the moments are corner shared 
between tetrahedra, the constraints lead to a correlated spin 
structure. The ‘disordered’ but long-range correlated state 
generates characteristic features in the magnetic structure fac¬ 
tor, appearing as “bow ti es’^^ or “pinch-points’ES indicative 
of power-law correlation^^®. 

In the semiclassical, ^ ^ 1, case the classical degeneracy 
is partially lifted by the zero-point energy of quantum fluctua¬ 
tions at harmonic order, but there remains an infinite manifold 
of degenerate collinear ground stated. Further quantum fluc¬ 
tuations at anharmonic order break the degeneracy between 
the various harmonic ground states, yet they leave out a mas¬ 
sive but nonextensive degeneracy (smaller than the harmonic 
ground state j^. 

The ground state for S = 1/2 is argued to be a quantum 
spin liquicffS®. This spin-singlet ground state has a finite 
energy gap for triplet excitations. The spin-spin correlation 
function decays exponentially with distance with a correlation 
length shorter than the lattice spacinj^. 


B. Pyrochlore Hubbard model 

Heisenberg interactions beyond nearest neighbor, induce 
transitions to various ordered phases (e.g., collinear, ne¬ 
matic and multiple-q orderj^SlT Also, easy axis anisotropy, 
long-r^e dipolar interaction, etc., lead to multiple-q ordered 
phase®. It has been argued that beyond the Heisenberg limit 
the half-filled Hubbard model on the pyrochlore lattice can 
be expressed as a highly frustrated intratetrahedral spin model 
with weak intertetrahedral perturbation®. This model has an 
exactly solvable Klein point, about which the ground state is 
a three-dimensional quantum spin liquid over an extended pa¬ 
rameter region. This spin-liquid state hosts massive spinon 
excitations, which are deconflned and move in all three di¬ 
mensions within the latticd^. For the Hubbard model on the 
pyrochlore lattice, the only work on the Mott transition, that 
we know of, suggests a transition from a semimetal to a spin 
liquid Mott insulato®. However, detailed properties near the 
IMT are not available. 


In the pyrochlore iridates R 2 lr 207 , both the R (rare-earth 
or Y) and the Ir live in two interpenetrating pyrochlore struc¬ 
ture. Their physics is mainly dictated by the 5d electrons of 
Ir, which have strong spin-orbit coupling and moderate Hub¬ 
bard repulsion U (due to large spatial extent of 5d orbitals). 
This strong spin-orbit coupling lifts the orbital degeneracy of 
5d electrons, reduces the bandwidth and leads to an effective 
single band description in terms of pseudo-spin jef / = 1/2 
states. In the absence of U, the ground state is a semimetal or 
topological insulator depending on the ratio of spin-orbit cou¬ 
pling and hoppinj^ (f). This picture remains unchanged for 
weak U/t. For strong U/t, Hartree-Fock calculation shows 
the system becoming an “all-in-all-out” magnetic insulator. 
Near the magnetic transition, a topological Weyl semimetal 
phase shows up®. A more elaborate cluster dynamical mean- 
field theory (C-DMFT) calculation confirms this scenaric®. 

In the pyrochlore molybdates R 2 M 02 O 7 , the physics is 
mainly governed by the 4d electrons of Mo. The presence 
of a trigonal crystal held and strong Hund’s coupling leads to 
band narrowinj^. Electron correlation then acts on a back¬ 
ground involving double exchange ferromagnetism compet¬ 
ing with superexchange antiferromagnetism. A Hartree-Fock 
calculation shows the possible ^in and orbital order in the 
ground state for the molybdate®. 


III. MODEL AND METHOD 

We study the single band Hubbard model, with nearest 
neighbor hopping, on the pyrochlore lattice: 

H = Ho + U^rii^ni^ (1) 

i 

where Hq = “ l^^ij)4cr<^Jcr- The Uj = -t for 

nearest neighbor hopping on the pyrochlore lattice and U > 0 
is the Hubbard repulsion. We will set t — 1. The chemical 
potential /r is varied to maintain the density at n = 1 as the 
interaction and temperature T are varied. 

We use a Hubbard-Stratonovich (HS) transformatior®that 
introduces a vector held mi(r) and a scalar held (tiji r) at 
each site to decouple the interaction. This decompositioi®® 
retains the rotation invariance of the Hubbard model, and 
hence the correct low energy excitations, and reproduces un¬ 
restricted Hartree-Fock theory at T = 0. 

We treat the m, and (pi as classical fields, i.e, neglect their 
time dependence, but completely retain the thermal fluctu¬ 
ations in m;. (pi is treated at the saddle point level, i.e, 
(pi — >■ {(pi) = {U /2){{ni)) = U /2 at half-filling, since charge 
fluctuations would be penalised at temperatures T <^U. Re¬ 
taining the spatial fluctuations of 111 ^ allows us to estimate 
scales, and access the crucial thermal effects on transport. In 
the literature the overall scheme is known as the ‘static path 
approximati on’ (S PA) to the functional integral for the parti¬ 
tion functior®131 jj^yg ^gg(J jjj ^jjg address the 

Mott transition on the triangula® lattice. Others have used 












it successfully in superconductorJ^S etc. We will discuss the 
limitations of this scheme later in the paper. 

Within this approach the half-filled Hubbard problem is 
mapped on to electrons coupled to the field m.i, which itself 
follows a distribution function 

U U 2 

Heff = ^ 2^ m, 

i i 

P{m,;} oc (2) 

where the chemical potential p in iJo is replaced by /i = 
/i — i?e// can be seen as comprising of an electronic 

Hamiltonian, Hei (the first two terms) and the classical ‘stiff¬ 
ness’ Hci = ^ ™ exact calculation, where the 

dynamics of riii and (j)i are retained, Hef / would be replaced 
by an effective action while P would be replaced by a fermion 
determinant in the {m, background. 

Within SPA i^e// and P{nii} define a coupled fermion- 
local moment problem. If the moments are large and ran¬ 
dom the electronic problem requires numerical diagonalisa- 
tion. Similarly, the P{m,} cannot be written down in closed 
form since the fermion free-energy is not known for arbitrary 
{iiii} background. The method of choice in these situations 
is a combination of Monte Carlo (MC) for updating the rrii 
with exact diagonalisation (ED) of the fermion Hamiltonian 
for computing the Metropolis update cost. 

To access large sizes within limited time, we use a cluster 
algorithm for estimating the update cost. The energy cost of 
updating the variable mi is computed by diagonalizing a clus¬ 
ter (of size Nc, say) constructed around the site R^. We have 
extensively benchmarked this method^. Results in this paper 
are obtained for lattices of 6^ unit cells with 4 atoms per unit 
cell, using a cluster with 3^ x 4 atoms. 

Electronic properties are calculated by diagonalising Hei 
on the full lattice for equilibrium {nii} configurations. Erom 
the equilibrium configurations we compute the single-particle 
density of states (DOS) as N (w) = ^ {J2n ^n)} where 
(...) indicates thermal average over equilibrium configura¬ 
tions. Eor a L X L X L pyrochlore lattice, the total number of 
lattice sites N = AL^. We calculate the optical conductivity 
by using the Kubo formuld^as follows, 

„XX(, _ ^0 / i^n) ~ /(c I rnm\2x(, , 

tr [U)) — / , _ \Jx \ \^m £«)) 

n,m 

where J""* = {n\Jx\m) and Jx is the current operator, given 

by, 

Jx — ^ ^ ^ H.C.) 

2,<7 

/(e„) is the Eermi function, e„ and |n) are the single¬ 
particle eigenvalues and eigenstates of P[e.i{mi] respectively. 
The conductivity is in units of ctq = /{Tiao), where 
oo is the lattice constant. The d.c. conductivity is ob¬ 
tained as a low frequency average of the optical conductivity 
over a window ~ O.OSf. We calculate average local mag¬ 
netic moment mavg = '^beir distribution 
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FIG. 1. Colour online: (a) Tight binding density of states for the 
pyrochlore lattice, (b)-(d) Results at T = 0 for: (b) variation of the 
average local moment mavg, (c) the density of states, A'^(O), at the 
Fermi level, and (d) the resistivity p{T = 0). 


P{m) = — |mi|)) from the equilibrium con¬ 

figurations. The magnetic structure factor is computed as 
‘S'(q) = at each temperature. The 

rapid growth of a mode at some q = Q (say) indicates the 
onset of magnetic order. Our results are averaged over 100 
equilibrium Monte Carlo configurations. 


IV. RESULTS 
A. The ground state 

As T —0 our MC mainly samples configurations that 
maximise P{mi}, or, alternately, minimize the energy, i.e., 
g^{Heff) = 0, This is the same as unrestricted Hartree- 
Fock in the magnetic channel. 

Upto a critical coupling Ud ^ 3.5f the minimisation yields 
rrii — I mil = 0 at all sites. As a result upto Ud the elec¬ 
tronic ground state is essentially tight binding, the density of 
states for which is shown in Fig[2(a). There is a sharply sup¬ 
pressed DOS at Fermi level characteristic of the pyrochlore 
band structure, and a flat band right above the Fermi level. 
For U < Ud the system is a semimetal. 

For Ud < U < UpG where Upo ~ 5t we observe a small 
moment, orientationally disordered, magnetic state. The aver¬ 
age moment size is shown in Fig[T](b), and the full distribution 
later in the paper. With the disorder caused by these moments 
breaking the translation invariance of the pyrochlore lattice, 
the DOS at the Fermi level, N{0) = N{uj)duj/ duj, 
where H = 0.05f, gains weight (Fig|T](c)). In a narrow re¬ 
gion around Upc there is rapid increase in the mean magni¬ 
tude of the moments and, as a result, the DOS at Fermi level 
gets depressed again. The detailed behavior of the low en¬ 
ergy DOS is shown later. At Uc 2 ~ 5.2f the DOS at the 












FIG. 2. Colour online: Phase diagram (left) and the density of states 
at the Fermi level (right) for varying U and T. Local moments appear 
at Uci but the state remains metallic, turning insulating at Uc 2 - PG 
refers to a pseudogap state and the metal-insulator transition line sep¬ 
arates regions with opposite signs of dp/dT. The right panel high¬ 
lights the ‘re-entrant’ feature in the low energy DOS with increasing 
temperature. To avoid clutter we have not marked the Upo scale in 
the ground state in the left figure. 


Fermi level vanishes as a Mott gap opens. For U > Uc 2 
the moments are large, and saturate to their atomic value, 

I mi I = 1, as(7/t—>oo. The coupling of the electrons to 
the local moments leads to weak scattering and a small resis¬ 
tivity for Uci < U < Upo, a rapid growth in resistivity for 
UpG < U < Uc 2 , and zero d.c conductivity for U > Uc 2 - 
This is shown in Fig[2(d). We will discuss the resistivity in 
much greater detail later, and just wanted to highlight the ef¬ 
fect within the T = 0 mean-field state here. 

A comment about the magnetic state. Since there was 
no reason to expect that the moments would have any obvi¬ 
ous periodic pattern the only way to do the ‘minimisation’ 
was via simulated annealing employing Monte Carlo. In the 
U/t ^ oo limit the half filled Hubbard model leads to a 
Heisenberg model for the (classical, in our SPA approx¬ 
imation) and for the pyrochlore Heise nberg antiferromagnet 
the mom ents a re known to be disorderecPSEIl albeit power law 
correlateci^lllll Our MC minimization reproduces this state. 
At lower U/t, both in the Mott phase and the ‘metal’, there 
are no simple results known - but our results suggest that a 
‘disordered’ state persists and correlations reminiscent of a 
classical spin liquid phasd^ survive down to U/t ^ 8. 


B. Thermal phase diagram 

Figure (left panel) shows the U — T phase diagram in 
terms of the magnetic, transport, and spectral properties that 
we observe. The following features emerge: 

At finite T thermal fluctuations of local moments on the 
weak-coupling side (U < Upc) lead to a quick low T in¬ 
crease in the low energy DOS, and then a gradual decrease 
with further increase in T. On the strong-coupling side (U > 
Uc 2 ) the angular fluctuations of the local moments result in a 
slight smearing of the Mott-gap with temperature and an in¬ 
crease in the low energy DOS. However in the Mott-transition 
neighborhood, the Mott gap quickly converts to a PG with in¬ 


creasing T, leading to the widening of the PG region shown 
in Fig. |^(left panel). 

We demarcate the finite T metal-insulator boundary in 
terms of the temperature derivative dp/dT. A state is metal¬ 
lic if dp/dT > 0 and insulating if dp/dT < 0. The spectral 
features and resistivity are discussed in detail further on. 

Figure (right panel) shows the DOS at the Fermi level 
varying with U and T. On the Mott insulating side {U > 
Uc 2 ) we observe the DOS slowly increasing with temperature, 
which can be understood as the filling of the Mott gap. On 
the metallic side (Ud < U < Uc 2 ) we see a non monotonic 
behavior. The DOS quickly grows with temperature in the 
low temperature regime. It then reduces with further increase 
in temperature, the weight getting transferred to high energy. 


C. Density of states 

Figure shows the thermal evolution of the DOS in three 
of the four broad interaction regimes of our phase-diagram. 

ii) For U < Uci the ground state is characterized by rrii = 
I mil = 0. The electron model reduces to the usual tight- 
binding pyrochlore lattice. This is characterised by two Hat¬ 
bands at the upper band edge and vanishing DOS at the Fermi 
energy. At finite T, small ‘randomly’ oriented local moments 
appear in the system broadening the flat bands and leading to 
a small DOS at the Fermi level, (ii) For Ud < U < Upo, the 
ground state has small disordered local moments. The DOS 
is gapless and the weight at the Fermi level is nonmonotonic 
with T, increasing initially and then decreasing as a weak PG 
forms (see Fig. |^(a)). {Hi) For Upc < U < U^ the DOS 
has a PG at T = 0. This fills up initially with increasing T, 
Fig. [^(b), but deepens again above a temperature scale that is 
visible in Fig. [fright panel, {iv) For U ^ Uc 2 the ground 
state has a hard gap. With increase in temperature, the angular 
fluctuations of the local moments result in a slight smearing 
of the Mott-gap, and an increase in the low energy DOS Fig. 
|3](c ). However, there exists a clear Mott-gap until very high 
temperature, T ^ niavgU. 


D. Transport and optics 

Figure 1^ shows the d.c. resistivity p{T) for different U/t. 
We consider the four regimes. 

(i) For U < Ud the T — 0 phase is a semimetal. Since 
this lies well below the Mott transition we do not show the 
T dependence here, {ii) For Ud < U < Upc the residual 
resistivity p(0) is finite with dp/dT > 0 over the entire T 
range. The resistivity can be understood in terms of a disorder 
induced density of states and the scattering of electrons from 
the small disordered moments. This is the metallic regime. 
{Hi) For U ^ Uc 2 the system has a clear Mott gap at T = 0 
with p(0) —>■ oo. In this regime dp/dT < 0 over the en¬ 
tire temperature window. This is the Mott-insulating regime. 
{iv) In the neighborhood of Uc 2 , i-e, \U — Uc 2 \ ^ Uc 2 , p{T) 
shows a non-monotonic behavior. We observe dp/dT < 0 in 
the low temperature limit, crossing over to dp/dT > 0 with 
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FIG. 3. Colour online: Density of states varying with temperature at three representative regimes of our calculation, (a) U — 0.85?7c2 lies in 
the gapless metallic side, (h) U = 0.9817^2 corresponds to the pseudo-gap (PG) regime with a pronounced dip in the DOS at the Fermi-level, 
and (c) U — 1.5Uc2 lies in the gapped Mott-insulating side. Tiow = O.Olf in panel (a), where as Tiow = 0.0 for panels (b) and (c). In panel (a) 
the low energy DOS reduces with increasing T, in panel (b) it increases and then reduces with T, and in panel (c) it monotonically increases 
with T. 


increasing T. The temperature at which dp/dT changes its 
sign is indicated as the Tmit- 

We observe Tmit increasing with U as seen in Fig|^ This 
behavior can be understood as the scattering of electrons from 
the background fluctuating local moments. As U increases, 
the average local moment magnitude mavg{U) also increases, 
resulting in the increased scattering of the electrons and a de¬ 
pleting DOS at the Fermi level. 

Figure shows the optical conductivity from our calcula¬ 
tion as we cross the Mott transition. The important points are 
as follows: (i) a{uj) for U < Uci is a semimetal at T = 0 
and does not have a Drude peak, (ii) For Uci < U < Upc, 
cr(uj) shows a response with the peak at a small finite fre- 
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FIG. 4. Colour online: Temperature dependence of the resistivity for 
different U/t near the Mott insulator-metal transition. The normal¬ 
izing scale is po = hje}. 



FIG. 5. Colour online: Optical conductivity at U/Uc 2 =0.85, 0.98 
and 1.5 with varying temperatures. U = 0.85C/c2 shows a non-Drude 
like behavior with peak at small and finite frequency. With increasing 
U the peak moves to higher frequency and the zero frequency weight 
decreases continuously and eventually there appears a gap for U > 

Uc2. 


quency that slowly shifts to higher values with increasing T. 

(iii) For U > Uc 2 the system has a clear gap A(T) in the 

DOS. Thus a{u}) = 0 for u) < ujc With increas¬ 

ing temperature the gap A{T) reduces, resulting in small, but 
increasing low frequency weight of a{uj) and the peak posi¬ 
tion shifts to higher frequency. This Mott-insulating regime 
of the pyrochlore lattice may have finite spectral weight at 
w = 0 in the optical conductivity a{uj) at high temperatures. 

(iv) For Upc < U < Uc 2 we have a pseudogap in the DOS. 
a{u: = 0) —^ 0 in this regime. However with increasing 
temperature the zero frequency weight increases initially and 
then decreases in accord with the behavior of the DOS in this 
regime. 


























V. DISCUSSION 

A. Overall scenario 

Within our picture, the interaction effects are encoded in 
the presence of the ‘local moments’ nii. The size rrii of this 
moment dictates the on site splitting at the site R^, and leads 
to a Mott gap in the overall DOS when Uruavg t. The 
spatial correlations of rrii decide whether the electron spin (7^ 
will display any long range order. 

Let us correlate the electron physics across the Mott transi¬ 
tion to the behavior of the nii, we will then take up the effec¬ 
tive models for the m, themselves, (i) In the metallic regime, 
U ^ [/cl, the mi’s are small and the orientations are random. 
The scattering from these moments leads to a finite broaden¬ 
ing, of the momentum eigenstates. This generates a finite 
DOS at w = 0, a non-Drude optical response, and resistivity 
increasing with temperature (unlike in a semimetal), (ii) In the 
Mott-insulating regime, U > Uc 2 , the mi’s are large and show 
short-range magnetic correlation (discussed next). The large 
mi’s lead to a Mott gap in the DOS and optical conductivity, 
and a diverging resistivity as T —>^ 0. (iii) In the pseudogap 
regime, the m^’s are moderately large and orientationally dis¬ 
ordered. This results in a strong suppression in the DOS at 
Fermi energy, but no gap, a large finite residual resistivity and 
a non-Drude optical response. 

We have shown niavg for varying [/ at T = 0 (Fig. [^b)). 
In the next section we discuss the limiting models that dictate 
the behavior of mi, and in the section after, we show detailed 
results on the size distribution of m^, and its spatial correla¬ 
tions. 


B. Effective models in limiting cases 

We have shown only a formal expression for the effective 
magnetic Hamiltonian. To get a feel for the magnetic states 
that arise it is useful to provide the approximate analytic struc¬ 
ture in limiting cases. These are (a) weak coupling, when 
U ^ t, and (b) strong coupling, when U ^ t. 

1. Weak coupling 

Our effective interaction looks like a ‘Hund’s coupling’, 
with the electron spin coupled to a background moment 
through the coupling U. Given the formal similarity with the 
Hund’s problem we can borrow the form of the weak coupling 
resulP21from the literature; 

. 1 

iTe//{m,}- - )m,.mj 

JjA 

+ 36 X! mfc, m3 -f .. 

ij 

The structure of /(m^, ...m;) is complicated, involving two- 
spin terms such as m^.mj, (m^.m^)^, three-spin terms 


such as (mi.mj)(mi.mfc), and four-spin terms such as 
(mi.mj)(mfc.m3. The i, j, k, I can be separated by long dis¬ 
tance in the weak-coupling limit. 

In contrast to the Hund’s problem, where there were pre¬ 
defined local moments, in our case the moments have to arise 
from an instability in the electron system. The leading in¬ 
stability involves the vanishing of coefficient of the quadratic 
term, i.e., 1 — [/xo(q) = 0. This generates the moment and, if 
there is a prominent peak at some q = Q in xo the moments 
order with that wavevector. In such a situation the fourth order 
term can be expanded about Q. The quartic term decides the 
magnitude of the order at Q. 

We have computed the Xo(q) for the half-filled pyrochlore 
band, having tested the scheme for the quarter-filled band for 
which results are availabld^. The Xo(q) is featureless, sug¬ 
gesting that there is no particular wavevector that would be 
picked out. In that case the quartic term, whose detailed struc¬ 
ture we do not know at half filling, decides not only the magni¬ 
tude but also the spatial character of the order parameter field. 
It seems that the nonlinearity creates a bimodal distribution 
for the rrii, discussed further on, but without any significant 
spatial correlation. 

2. Strong coupling 

For U ^ t, one can write an effective magnetic Hamilto¬ 
nian on the pyrochlore lattice by tracing out fermions order by 
order in t/U. This gives us 

HtetrittO-} ~ Jq + JlPa + J + •■•) 

Oi 

fcGan/3 

iFgg„p{m}~4 ^ S,.S,+4' ^ (S,.Sfe)(S,.Sfc)... 

i£o6,jG^ i£cx,jG0 

where spin on the tetrahedron 

a, Htetr describes interactions between spins in a tetrahe¬ 
dron while Hcoup includes the intertetrahedron terms with a 

common corner shared site. Jq = —^{1 — J 2 = 

1,(1-^), J4 = i;andJ',4'^0(^). 

Deep in the Mott phase, one would drop the J 4 , J4 and ,]'( 
terms and obtain a classical Heisenberg model, which does not 
show any long-range or der or freezing, but power-law correla¬ 
tions at low temperaturdSEI. The correlations of the electron 
spins, (Ji, can be computed in response to this. As U ^ Uc 2 
the J 4 , J 4 and J'l terms become important. These multi-spin 
exchange interactions modify the magnetic ground state, but 
some of the Heisenberg limit featured are observable (see 
next) in the structure factor down to [/ ~ 8f. The expansion 
in t/U ceases to be useful once the gap closes. 

C. Detailed magnetic structure 

Figure]^ shows the the amplitude distribution P{m) of the 
magnetic moment for different T and interaction regimes. 





FIG. 6. Colour online: Temperature dependence of P{m) for U = 0.85?7c2, Uc 2 and 1.5?7c2 for indicated temperatures. 
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FIG. 7. Colour online: The full magnetic structure factor 5(q) for U = 0.85(7c2, Uc 2 , 1.5?7c2 (along column) and T = 0.0, 0.02f and O.OSf 
(along row). We use the notation q = ^ {ux, Uy, n^), n’s are integers. The size of a dot signifies relative weight at a given q while its color 
represents the actual magnitude of S'(q). 


Figure|^a) shows Pirn) for U = 0.85(7c2. In the ground 
state there is a two peak structure, highlighting the presence 
of amplitude inhomogeneity. We have checked that there is no 


signihcant density inhomogeneity, or charge ordering, in the 
system. With rise in temperature, P{m) shows a broad single 
peak behavior with the peak shifting towards large m. This 



























































behavior is seen in the Uci < U < Uc 2 regime. Figure l^b) 
shows P{m) for U ^ Uc 2 , just in the Mott insulating side. 
The ground state has a narrow single peak feature indicating 
an amplitude homogeneous Mott state. With rise in temper¬ 
ature this narrow peak broadens, the peak position moves to¬ 
wards higher m. Figure l^c) shows P{m) for U = 1.5Uc2, 
well in the Mott regime. P[m) has a single peak feature which 
broadens with temperature and shifts to higher m. The fluctu¬ 
ations about the mean are weaker in the insulator than in the 
metallic state. 

Figure shows the q dependence of the magnetic struc¬ 
ture factor for varying T and U. We observe that even at 
T = 0, 5'(q) has no ordering peak at any q’s. The magnetic 
ground state is disordered. Nevertheless for both U ~ Uc 2 
and U = 1.5C/c2 the weight distribution is not completely 
homogeneous in q and have some prominent features. This 
signature survives to T ~ 0.03f. Examination of the classi¬ 
cal Heisenberg model has revealed that there are ‘pinch point’ 
features in scattering, arising from the constrained ground 
state, that survivd^to T ~ O.IJ, where J is the exchange 
scale (in our case J = /U). Obtaining such a result, in¬ 

dicative of power law correlations, requires larger system size 
(> 10^), much longer annealing (> 10®) MC weeps), and 
a more sophisticated algorithm instead of single spin update. 
Due to our computational cost and size limitations we only get 
a hint of this spin liquid state. By the time T ~ O.OSf, which 
is ~ J/2 at 17 = 1.5 C/c2 , the q dependence is completely 
featureless. 


D. Methodological issues 

There are several approximations, analytic and numerical, 
we had to make to achieve some headway. Let us comment 
on these one by one, and their effect on the reliability of our 
results. 


1. Static approximation in the Monte Carlo 

In principle both spatial and temporal fluctuations could be 
important near the Mott transition on the pyrochlore lattice. 
However, fully handling the temporal fluctuations of the mi 
and ())i fields requires a quantum Monte Carlo scheme. The 
absence of such results is probably due to the sign problem 
for fermions in the frustrated geometry. We have retained only 
the zero Matsubara frequency, = 0, mode of the auxiliary 
fields, exactly. 

It is possible, but non trivial, to set up a Gaussian expansion 
for the finite modes of and (jti, while retaining and 
treating the = 0 mode exactly as we have. This is a project 
for the future. In the absence of such a calculation we can only 
make the following conjectures. 

(i) Local moments would remain well defined in the insu¬ 
lating phase, due to the gap. Since the present theory itself 
predicts no order at any U or temperature, we do not ex¬ 
pect additional quantum fluctuations to qualitatively modify 
the magnetic state, (ii) A possible qualitative consequence 



FIG. 8. Colour online: Comparison of the variationally obtained 
average moment value with result from the Monte Carlo as T —>■ 0. 

of quantum fluctuations could be the restoration of translation 
invariance in the ‘magnetic metal’ phase. We cannot make 
a definite comment on this but even if this were to happen, 
making the metallic state perfectly conducting at T = 0, we 
believe that above a low coherence temperature one would see 
the signature of a highly resistive metal. 

2. Finite size effects 

We have done our calculations on lattices that involve ^ 
800 atoms. This is much larger than what is typically ac¬ 
cessible in full-fledged fermion quantum Monte Carlo (QMC) 
(since temporal fluctuations also have to be included), allow¬ 
ing us to access resistivity and spectral features which do not 
involve artificial broadening. It is possible to try the calcu¬ 
lations on 8® X 4 lattices (involving about 2000 atoms) and 
checks at individual parameters do not reveal a significant dif¬ 
ference. 


3. Cluster based update 

Our Metropolis update involves a small cluster rather than 
diagonalisation of the full Hamiltonian. This is well con¬ 
trolled in the large U limit when the ‘range’ of electron ex¬ 
cursion is limited but less reliable near the Mott transition. 
Lor that purpose we have used a variational calculation, de¬ 
scribed below, that uses the MC result as an ansatz and checks 
the stability of such a state on large (up to 10® x 4 lattices) 
for the ground state. The results are discussed below - and are 
qualitatively consistent with the cluster based MC. 

4. Variational check 

Since MC hints that the magnetic ground state involves 
disordered local moments for U > Ud, we tried a simple 
variational check. We set up trial configurations {riii} with 
random orientation, but uniform magnitude niQ, with rriQ as 
a variational parameter. This differs from the real situation 
where the mi’s have some amplitude inhomogeneity and also 



Indicator 

Present theory 

Molybdate 

Iridate 

Degrees of freedom 

One orbital per site 

Two orbitals per site and a. S — 1/2 
moment 

Effective one orbital per site 

Interactions 

Repulsive Hubbard 

Repulsive Hubbard, Hund’s coupling, 
antiferro superexchange 

Repulsive Hubbard, spin-orbit 
coupling 

Weak-coupling phase 

Band semimetal, spin disordered 
metal 

Ferromagnetic metal, spin-glass metal 

Correlated metal 

Strong-coupling phase 

Spin-liquid Mott insulator 

Spin-glass Mott insulator 

All-in-all-out Mott insulator 

Features in resistivity 

Continuous divergence of p(0) with 
Vjt 

Seemingly continuous divergence of 
p{T = 0) with rn 

Seemingly continuous divergence 
of p(T = 0) with rn 

Magnetoresistance in mod¬ 
erate field 

Weak 

Huge: drives insulator-metal transition 
(Gd2Mo207) 

Weak 

Anomalous Hall response 

Not explored 

Observed in Nd 2 Mo 207 

Observed in Nd 2 lr 207 


TABLE I. Comparison of features of the pyrochlore Hubbard model with the models relevant for the molybdates and iridates. 


orientational correlation. Energy minimisation confirms the 
presence of a small moment phase, beyond some Uiow, with 
an initial slow growth of mo(C/) and then a rapid crossover 


E. Relation to experiments 

Pyrochlore Mott materials include the rare-earth based 
molybdates R 2 M 02 O 7 and iridates R 2 lr 207 (R is a rare-earth 
ion or Y). The Mo and Ir ions live on a pyrochlore lattice while 
the R inhabit an interpenetrating pyrochlore structure. 

In contrast to 3d electron based systems which are domi¬ 
nated by the Hubbard U, the molybdates involve 4d electrons 
- where the Hund’s coupling Jh is also important, while the 
iridates involve 5d electrons - with spin-orbit coupling playing 
a vital role. In this paper we focused on the simple Hubbard 
model as the starting problem. In what follows we quickly list 
out the key molybdate and iridate Mott features (see Table |I]i 
and argue how some of these are already visible in our present 
results. 

The molybdates exhibit a MIT with decreasing rare-earth 
ionic radius ta- Materials with larger va are ferromag¬ 
netic (F M) me tals, those with small ta are ™in glass (SG) 
insulator^i^ini and there is a SG metal phas JSI near the MIT. 
Pressure and magnetic field can drive an IMT in materials 
like Gd 2 Mo 207 which are weakly insulating^. An anoma¬ 
lous Hall effect (AHE) has been observed^^l^ in Nd 2 Mo 207 
and is ascribed to non vanishing spin chirality. 

The iridates also show a MIT with decreasing ta, but in 
this case the MI transition is accompanied by a magnetic tran¬ 
sition from a paramag netic to an antiferromagnetic ‘all-in-all- 
out’ (AIAO) orderin^i^ISIlII While the magnetic character 
differs distinctly from molybdates, iridates also show a pres¬ 
sure driven i nsulato r-metal transitions, via unusually resistive 
ground state ^pin chirality driven AH^^in materi¬ 

als like Pr 2 lr 20 ,ISl. 

Reducing va reduces the hopping in these materials and 
is akin to incre asing the U/t ratio in our model. The gen¬ 
eral observatiorP^EEIlIll of a metal-insulator transition with 
reducing ta is consistent with our phase diagram. In the 


to large values at some Uhigh (see Figure]^. The Uiow and 
Uhigh are about 10% higher than MC estimates and may get 
reduced if spatial correlations are included. 


table above, however we compare a set of indicators be¬ 
tween molybdate or iridate experiments and our present the¬ 
ory. While the presence of electron correlation on a py¬ 
rochlore structure is common to all three, there are additional 
interactions present in the real materials. The one common 
feature that emerges is the occurence of a high resistivity 
ground state close to the MIT which, we believe, is due to the 
coupling of electrons to disordered spin and orbital moments. 
This is a general consequence of the interplay of correlation 
effects and strong geometric frustration. We will discuss the 
detailed iridate and molybdate theories elsewhere. 


VI. CONCLUSION 

We have studied the Mott transition in half-filled Hubbard 
model on a pyrochlore lattice. The geometric frustration and 
the corresponding large magnetic degeneracy prevents the oc¬ 
curence of any magnetic order in the deep Mott state. This 
continues all the way to the insulator-metal transition. Be¬ 
yond the insulator-metal transition there is a window with a 
pseudogap in the density of states, disordered local moments, 
and a large residual resistivity. At even weaker interaction 
one recovers the non magnetic band semimetal. Thermal fluc¬ 
tuations destroy the ‘spin-liquid’ correlations in the insulat¬ 
ing state, converting the system to an uncorrelated paramag- 
net. The low energy electronic density of states and the re¬ 
sistivity show a monotonic temperature dependence deep in 
the metallic and insulating phases, but a non monotonic char¬ 
acter near the insulator-metal transition. A comparison with 
the rare-earth molybdates and iridates reveal a clear similarity 
in the resistivity, while detailed aspects of the magnetic state 
and magnetotransport are missed out. We are exploring these 
problems separately. 
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